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Abstract 

The  linearized  theory  oi  liee-auriacc  gravity  i lew  past  sub¬ 
merged  or  floating  bodies  is  based  on  a  perturbation  expansion  of 
the  velocity  potential  in  the  slenderness  parameter  f  ,  while  the 
Froude  number  V  is  kept  fixed.  It  is  shown  that  although  the 
free  waves  amplitude,  and  the  associated  wave  resistance,  tend  to 
zero  as  F  +  0  ,  the  linearized  solution  is  not  uniform  under  this 
limit:  the  ratio  between  the  second  order  and  first  order  terms 

becomes  unbounded  for  F -+  0  and  fixed  c  .  This  non  uni  tornu  t  y 
(called  "the  second  Froude  number  paradox"  in  a  previous  work)  is 
related  to  the  nonlinearity  of  the  r ree-surf ace  condition.  Cri¬ 
teria  of  uniformity  of  the  thin  body  expansion,  combining  t  and 
F  ,  are  derived  for  both  two-  and  three-dimensional  flows.  These 
criteria  depend  on  chi  shape  ot  the  leading  (and  trailing )  edge*: 
as  the  shape  becomes  finer  the  linearized  solution  becomes  valid 
for  smaner  F  . 

Uniform  first  order  approximations  are  derived  by  two  alter¬ 
native  methods:  velocity  straining  and  coordinate  straining. 

In  the  first  case  the  uniform  unperturbed  velocity  in  the  tree- 
surface  condition  is  replaced  by  a  variable  velocity  distr: but  ion . 
The  second  method  leads  to  an  apparent  displacement  of  the  most 
singular  points  of  the  body  skeleton.  In  both  cases  the  para¬ 
meter  v.  appears  not  only  in  the  wave  amplitude,  as  implied  by 
the  tli in  bouv  expansion,  but  also  in  the  wave  number  function. 

The  nonuniformity  of  the  usual  thin  body  expansion  is,  there  lore, 
similar  to  that  encountered  in  problems  character ized  by  multiple 
scales . 


. . . 
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INTRODUCTION 


The  linearized  theory  ot  free-  aurtacc  gravity  flow  past  sub¬ 
merged  or  loafing  bodies  is  based  on  the  assumption  that  the 
body  causes  a  small  disturbance  ol  a  uniform  flow.  Such  an  appro¬ 
ximation  is  incorporated  in  a  systematic  asymptotic  expansion  of 
the  velocity  potential  by  assuming  that  c  (beam  length  ratio 
for  a  thin  ship,  draft /length  ratio  for  a  flat  ship,  body  length 
submergence  depth  ratio  in  the  case  of  deep  submergence)  tends  to 
zero  while  the  Fronds  number  F  (hanofl  n;i  >'ody  length  o*"  sub¬ 
mergence  depth . respectively )  is  kept  fixed. 

In  previous  works  (Salvesen,  1969  ;  Dag  an,  1972a)  it  lias  been 
shown  that  it  is  not  legitimate  tc  let  F  -»  0 ,  for  a  fixed  e  ,  in 
the  linearized  solution,  or  in  oilier  words  that  the  usual  approxi¬ 
mation  is  not  uniform  in  F.  Two  "small  Froude  number  paradoxes" 
have  been  formulated  in  this  context  (Pagan,  1972b)  and  ad-hoc 
uniformizotion  procedures  have  been  suggested  (Ogilvie,  1968; 

Dagan,  1972b),  loading  to  a  quasi-linearization  of  the  frec-surface 
condition.  It  has  been  proved  (Tuck,  1965;  Salvesen,  1969; 

Pagan,  1972a)  that  the  small  Froude  number  nonuniformity  is 
associated  with  the  nonlinearity  of  the  free -surf ace  condition. 

In  all  cases  detailed  computations  have  been  carried  out  only 
for  two-dimensional  flows. 

In  the  present  study  the  problem  oi  the  small  Froude  number 
nonuniformity  is  attacked  in  a  different  way,  the  results  being 
alsr,  different  of  those  obtained  previously.  For  the  first  time 
the  influence  of  the  bluntness  of  the  bow  on  the  small  F  solution 
is  discussed  in  detail  and  the  analysis  is  extended  to  three- 
dimensional  flow  in  general  and  to  flow  past  thin  ships  in  parti¬ 
cular. 
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It  is  worthwhile  to  mention  here  that  the  problem  .is  related 
mainly  to  throe-dimensional  applications,  since  a  large  class  ol 
ships  operate  at  relatively  len.  1'roude  numbers  and  in  most  such 
cases  the  usual  theory  of  wave  resistance  has  been  found  to  be 
unsatisfactory.  We  begin,  nevertheless,  with  studying  the  two- 
dimensional  flow  because  the  use  of  the  powerful  tool  of  analyti¬ 
cal  functions  in  this  case  pennits  to  clarify  seine  matters  of 
principle  much  easier  than  in  three-dimensions . 

Obviously,  there  are  various  possible  factors  related  to  the 
discrepancy  between  wave  resistance  as  measured  in  experiments  and 
as  predicted  by  the  linearized  theory',  like  viscous  effects  or  the 
bow  breaking  wave.  This  should  not  deter  us,  however,  irom  seeking 
a  consistent  solution  for  the  wave  resistance  in  the  frame  ol  the 
potential  flow  theory. 
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PART  I 

TWO-DIMENSIONAL  FLOW  PAST  SUBMERGED  BODIES 


1 .  The  Thin  Body  Expan a ion 

We  consider  a  steady  uniform  flow  from  infinity  past  a 
submerged  body  (Fig.  1).  Let  z'  =x'  +  iy '  be  h  complex  variable, 
w’  =u’  -  i  v '  the  complex  velocity,  f  1  =  tp  '  +  ip'  the  complex  poten¬ 
tial,  n*  the  free  surface  elevation,  2L’  the  body  length,  h'  the 
submergence  deptn ,  2T'  the  thickness  ana  U’  the  velocity  of  uni¬ 
form  flow.  First,  variables  are  made  dimensionless  by  referring 
them  to  L '  and  U'  ,  i.e.,  z  ~  z' /!,'  ;  f  =  f!/L'  ;  w-w'/V  ; 
h=h'/L’  r  c  -  T '  /  L 1  a:  id  F  -  U  '  /  IgL’ )  1/2  . 

Under  an  expansion  of  the  analytical  function  f(z;e,h,F) 
in  a  small  o  asymptotic  series 


f  =  -  z  +  c  f  ^  ( z ;  h  ,  F )  +  i  ( z  ;  h  ,  F )  +  .  .  . 


(1.1) 


the  following  sets  of  equations  are  obtained  for  and 

from  the  expansion  of  the  exact  equations  (Wehausen  and  Laitone, 

1960) 


lm(iF‘ 


df 


<J: 


1 


n,  3  '!■ 


) 


(y  =  0 ) 


il.  z  i 


(]  .  3! 


h  *  U 


(x  -*•  00  ;  y  -*•  -m ) 


(1  .4) 


=  +t  (x) 


( j  x|  <  1  ,  y  =  “hi 0 )  (1.5) 
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4. 


where  n  --  er^  +  +  ...  and  y  =  -htct(x)  is  the  equation  of 

the  body  profile,  assumed  to  be  symmetrical  for  the  sake  of 
s implieity  ? 

df 


In  {iF‘  -r~-  f_)  =  (x)  =  -Fz  (i  u  *  +  i  v*)  +  F“u 


dz 


r  ui  7  vr 


?u 

.  4  ..  1 


‘1  Tx“ 


n-  =  ^_  +  t^u 
^  z  .1  ± 


f2  -  0 


(l.b) 

(y»0) 

(1.7) 

(1.8) 


(x  -r  °°  ;  y  +  -<*>) 

K'2  -  ?u3t  (|x|  <  1  ,  y  =  -htO)  (1.9) 

In  addition,  a  Kut ta- Jonkovsky  condition  has  to  be  imposed 
in  the  case  of  a  sha] 
unique . 


D  = 


trailing  edge  in  order  to  make  circulation 

ialvesen,  13G3)  that  far  behind  the 

;  the  expressions 

body 

Im  (a, 

X 

e  "x)  (x  -*■  -«) 

(.1.10) 

Im  (a., 

e  lx)  +  const  (x  ■+•  -») 

(1.11) 

given 

by  the  following  expression 

1  | 

2Y7  1 

2  1  2 

ca,  f  e*a2 I 

(1.12) 

!L‘  . 

Hence,  by  expanding  D 

eiDi 

+  c5D2  +  0\zk) 

(1.13) 

* 

'I 


-if 
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we  have 


Di  =  TfT 


D2  -  p-  Re  (^^2 ) 


(1.14) 


The  method  of  determining  f,  and  f7  ,  solutions  of 
Eqs.  (1.2)  -  (1.8)/  is  well  known.  Let  w^  and  w^  be  the 
first  order  lineari2ed  solutions  of  the  velocity  of  flow  past 
tne  body  or  its  image,  respectively,  in  an  infinite  domain, 
i  .e.  , 


w^(z)  = 

e 

l-ih 

/  t  (x  ) 

-l-ih  s 

dz 

s 

7T 

2  -  Z 

S 

where 

z 

s 

=  xs  + iys  is 

the 

coordinate 

of  a  point 

along  the 

skeleton 

('|xsl  <  1,  ys  = 

-h)  , 

Re  x  -  dt/dx  end 

( Z )  - 

i. 

l  +  ih 

/  T(X  ) 

-1+ih  8 

(z)  . 

tT 

W1 

2  ~  zs 

Then, 

the 

solution  for 

may  be  written  as 

(1.15) 


(1.16) 


where 


f,  =  f5  +  f?--  /  wV  (o)ui(~)cla  (Im  z  <  Im  o )  (1.17) 

1  J  1  71  •*  J.  f 


■  _  °°-io  i\ 

w(C)  »  e"1’  /  dA  =  / 


"  elp 


P  +  K 


dp 


(1.18) 


the  A  plane  being  cut  along  im  A  -  C ,  Re  A  >  0 ,  and  the  p 
plane  is  along  Jm{p+0  =0,  Re  (p+O  >0. 

The  second  order  solution  satisfying  (1.6)  and  regular 
in  the  lower  half  plane  (we  consider  here  only  the  f ree-surface 
second  order  effect  and  disregard  (1.9),  which  leads  to  less 
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severe  effects  as  F  -*  0)  may  be  written  as 


.  uo 

f  2  (2 )  =  (0 )  u)  (~rr)  d  0  (Imz  <  Im  a)  .  (1.19) 


2 .  The  Second  Order  Solution  (f ree-surf ace  effect) 

We  are  going  now  to  transform  (1,19)  such  that  f^  will  be 
expressed  as  an  integral  over  analytical  functions  of  0  .  First, 
we  have,  by  integration  by  parts 


r 

j 

•  to 


3v 


1  w(££)do  = 


'l  So-  'T-' 


(ulvi  +  vi* 


z  -  a 


(1.20) 


As  far  as  the  free  waves  are  concerned,  the  last  integral  in  (1.20) 

renders  the  well-known  Stokes  second  order  waves  of  amplitude 
-2h/F 2 

0(e  )  Under  the  limit  F  *»•  0  these  waves  are  negligible, 

/p2 

as  compared  to  the  remaining  terms  which  are  0(e  '  ),  and  will 

be  neglected  in  the  sequel. 

Hence,  by  (1.19)  and  (1.20)  we  have  under  these  conditions 


s  7  I  <; 


vl- 


3 

T 


U 


-iuivl) 


w(^r)dc 


(1.21)  . 


and  v^  .which  are  obtained  from  (1.17),  may  be  written  along 
the  real  axis  as  follows 
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Substituting  (1.22)  and  (1.23)  m  (1.21)  and  integrating  by 
parts  we  obtain  f2  (for  details  see  Appendix  I)  in  its  final 
form  as 

.  00 

f2{z)  ~  ?  /  *  “  *wl?(cr)  J2  + 

—  00  A  1  1 

0  00 

+  ~T  .fwJ(a)  +w^(o)  ]  /  w“(t)a)(^.)dx  }w(2l^)da  (1.25) 

Besides  the  term  containing  the  Stokes  second  order  waves  we  have 

neglected  in  deriving  (1*25)  also  the  term  related  to  the  constant 

part  of  (1.11),  which  is  associated  with  the  DC  part  of 

P-,  (x) ,  and  does  not  contribute  to  the  wave  res  i.  s  ti3jn.cc  • 

*» 


i 

1 

j 
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3.  Illustration  of  Results:  The  Source-Sink  Body 


( i )  Th t  oong  Zsta 

Rather  than  pursuing  a  general  discussion  of  the  second  order 
solution  under  the  limit  F-'O,  we  begin  with  a  simple  case  which 
can  be  solved  in  a  closed  form. 


We  consider  a  closed  body  genera  tec  by  a  source  at  53  1  -  ih 

and  p.  sink  at  z^_  =  -1  -  ih.  In  view  of  our  interest  in  three- 
dimensional  applications  we  consider  only  the  thickness  effect 
(the  influence  of  circulation,  generated  by  distributed  or  con¬ 
centrated  vorticity,  may  be  analyzed  in  a  similar  way) .  This  is 
the  first  order  representation  of  a  straight  thin  body  with  blunt 
leading  and  trailing  edges. 


Following  now  (it 15), 


(1.16),  (1.17)  and  (1.2b)  we  have 


f 


1 


z-z 

z-"z 


L 

t 


+ 


1 

TT 


1  Z"Zt 
7  “(-p T“> 


(1.26) 


i 

i 

s 


J 


i  _  ii  ii 

wi  -  j?  tt;  “  57  7^  ? 


(1.27) 


f 


2 


CO 


1 

7 


o-z, 


+ 


i 


i 

1 


x  O  ( ~pT§*)  dp }  w  (  3o 


(1.28)  a 

4 

4 
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The  integration  in  (1.28)  can  be  carried  out  exactly.  The 

t  u  i 

first  term,  resulting  from  (w^-w^)  ,  contributes  by  the 

residues  at  0  «  zx  and  o  *  z  .  The  last  terms  are  more 

i  t 

intricate,  but  still  tractable,  at  least  for  x  +  -°°.  ! 

We  consider  now  the  expansion  of  f^  (1.26)  and  f2  (1.28) 
for  small  F  . 

! 

(ii)  Th*  ~nsar  _£i0ld  _8oliition 

For  F2  -*-Q,  o)(z-z  /F2 )  can  be  expanded  in  an  asymptotic 

fL 

series  for  fixed  r  as  follows 


a)  ( 


y  (n-1) I  F2n 
n=l  in(z-z^)n 


(1.29) 


This  expansion  is  valid,  however,  only  for  Iz-z^l  >  6  and 
-  v  +  6  <  arg(z-z^)  <  0,  where  <5  and  3  are  arbitrarily  small 

fixea  quantities  (for  details  see  the  discussion  of  the  related 
exponential  integral  function  in  Copson,  1965). 

Substitution  of  (1.29)  and  the  similar  expansion  of 
[z-zt/F2]  into  w^  ,  obtaining  by  differentiation  of  f^  (1.26), 
yields 

=  w*  +  +  0(F?)  (-tt+2  <  arg(z-z^)  <  0)  (1.30) 

Hence,  w^  degenerates  at  zero  order  into  the  rigid  wall  solution, 
i.e.,  the  solution  of  flow  past  the  body  in  the  presence  of  a  rigid 
wall  at  y*0.  However,  this  limit  is  not  uniform  and  in  parti¬ 
cular  is  not  valid  far  behind  the  body,  i.e.,  for  x  -*•  -  00  and  y 


■lEwrunu 
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-2t/F2  -iz/F2  | 

supplemented  by  the  term  2trie  e  which  represents 

precisely  the  trailing  waves.  For  this  reason  (1.30)  may  be 

called  the  nearfield  expansion 

i 

The  rigid  wall  solution,  and  the  subsequent  terms  of  (1.  ’>0), 
may  be  obtained  also  by  expanding  first  the  linearized  free  sur¬ 
face  condition  (1.2)  for  F2  ■*  0  and  solvi’-'o  term  by  term.  In  j 

contrast  with  the  previous  procedure,  howevei ,  (1.30)  is  thus  ! 

obtained  as  a  uniformly  valid  solution  in  the  entire  z  plane.  ! 

This  difference  in  results,  manifest  in  the  lower  half  plane 

\ 

y  <  0  behind  the  body,  has  been  called  "the  first  small  Froude  j 

number  paradox"  in  a  previous  work  (Dagan,  1972b).  Although 

the  wave  term  is  exponentially  email  for  y<h,  as  compared  to 

the  powers  of  F2  in  (1.30),  it  is  the  only  one  which  does  not 

tend  to  zero  for  x -*•-«,  y  =  fixed  and  which  is  associated  with 

wave  resistance. 

Similarly,  the  near  field  expansion  of  f2\z)  may  be  ! 

obtained  frc.  (1.28)  for  F2  -*■  0  by  expanding  u)(z-a/F2)  and  ; 

computing  the  residue  contributions  at  and  zfc  .  The 

result  is  0(F2)  and  is  a  rational  function  of  z  with  poles 
of  different  orders  at  z  =  z^  and  z  =  z_  .  Hence,  ef^  is 
0(e)  while  e2f2  * 3  0(e2F2)  and  the  near  field  expansion  of  j 

f  is  uniform  as  F2  0 .  j 

i 

* 

I 

(lit)  The  £ree  w »s*s  jjOtvHti-uZ 

The  free  waves  potential  is  obtained  in  (1.26)  and  (1.28) 
by  .letting  x  ■*  -  00  .  The  firnt  order  solution  (1.26)  yields 


iz?/F2 

2i(e 


•4  (sin 


1  ,  -h/F2  -iz/F2 

ft)  e  e 


(1.31) 


»n),«an  Kummm*-  . . . 
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In  the  second  order  solution  (1.28)  w(z-o/F2)  is  first  replaced 

-iz/F*  io/F* 

by  2wi  e  '  e  '  Integration  yields  (for  details  see 

Appendix  II) 


w  1  — iz/F2  -h/F2 

f 2  -  "  fT  t(b2+ic2)e  *  +  (b2-ic2)e  *  ]e  1Z/F  +0(e  h/F  ) 


(1.32) 


where  b2  “  ~  (y  +  2C  +  In  4  +ln  »  c.*2  «  i  (y-arc  tg  i)  and 

C  is  Euler  constant.  If  h<<l,  c2  =  0  and  (1.32)  becomes 


(cos 


e"h/F2  e“iz/F2 


0(e"h/F2) 


(1.33) 


Hence/  the  amplitude  of  the  free  waves  has,  by  (1.10)/  (1.11), 
(1.31)  and  (1.32)  the  orders 


ea^  -  0(ee”h'/F  ) 
e2a2  =  0(e2e~h/F  /F2 ) 


(1.34) 

(1.35) 


and  although  for  F2  0  and  h  and  e  fixed  both  ea^  and 
e2a2  tend  to  zero,  their  ratio  becomes  unbounded  like  e/F2 . 

This  ncnuniformity  of  the  thin  body  expansion  has  been 
called  in  a  previous  work  "the  second  small  Froude  number  paradox" 
(Dagan,  1972b)  and  it  has  been  described  previously  by  Salvesen 
(1969).  Eq.  (1.35)  shows  that  the  usual  linearized  theory  is 
valid,  for  the  source-sink  body,  only  if  e/F2  =  o(l),  i.e., 
for  large  Froude  numbers  based  on  thickness. 
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4 .  Generalization  for  Bodies  of  Different  shapes 

0  ij 

Since  for  an  arbitrary  thickness  distribution  and 

are  repreaented  by  source  distributions  (1.15),  (1.16)  the  results 

of  the  previous  section  may  be  extended  to  thin  bodies  of  any 
shape.  It  is  easy  to  ascertain  that  the  near  field  solution, 
based  on  (1.17)  and  (1.29),  has  the  rigid  wall  approximation  as  a 
leading  term  and  is  uniform  in  the  sector  tt  ~  6  <arg(z-l-ih)  <  0 
as  F*  0. 

The  nonun if ormity  of  the  expansion  of  the  free  wave6  depends 
essentially  on  the  bluntness  of  the  leading  edge  (for  the  sake  of 
simplicity  we  consider  bodies  of  a  smooth  shape  and  assume  that 
viscous  effects  ensure  anyway  that  the  trailing  edge  has  a  fine 
shape).  The  free  waves,  at  first  order,  are  represented  by 

_w  _  -iz/F2  I  u,  >  ia/F2, 

»  2e  J  (a)  e  do  (1.36) 


which  ha6  been  obtained  from  (1.17),  the  integration  path 
circumventing  the  skeleton  of  the  image  of  the  body  -1  <  a  ~ih  <  1 
in  the  upper  half  plane.  For  F2  -*•  0  the  integral  in  (1.36) 
may  be  expanded  in  the  usual  manner,  the  lowest  order  term  bei«:g 
provided  by  uhe  highest  singularity  of  w^(o),  at  a  =  1  +  ih. 


ax  =  0(e'h/F‘)  . 


We  have  seen  that  for  a  source -like  blunt  shape 
For  an  elliptical  shape,  (i.e.  wY  v  l//a-l-ih  )  of  the  leadi 

r  /■—:  v 

edge,  (1.36)  shows  that  a.  =  0(Fe*"n/r  ).  Similarly,  for  a  wedge 


na 


u 


like  shape  (w 
(see  Lighthill,  1964). 


^  'v  ln(o-l-ih))  we  obtain  a^-0(F  In  Fe 


^"h/F‘v 
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To  estimate  the  order  of  the  amplitude  of  the  free  waves 
at  second  order  we  have  to  use  the  expression  of  f_  (1.25), 
with  w(z>o/r  )  replaced  by  27Tie  t'/ L  o  .  The  computation 

is  facilitated  by  the  observation,  supported  by  the  detailed 
solution  of  the  previous  section,  that  the  order  of  the  lowest 
term  in  F  is  determined  by  the  term  [w^(c)]2  in  the  integral 
of  f2(z)  (1.25),  the  other  terms  contributing  at  an  equal  or 

higher  order.  Hence,  the  order  of  is  determined  by  integrals 

of  the  type 

e'iz/F*  j  Iwu(o)J2  d0  .  (1.37) 

We  have,  therefore,  for  an  elliptical  leading  edge 
a2  -0(e  h/^  )  and  for  a  wedge-like  shape  a2  =  0(F2ln2F  e  il/ F  ). 
In  each  case  the  far  waves  amplitude ,  and  consequently  the  wave 
resistance,  is  not  uniform  for  F2  -*-0,  the  nonuniformity  becoming, 
however,  weaker,  as  the  shape  of  the  edge  becomes  finer.  The 
results  are  collected  in  the  following  table: 

TABLE  1 


The  shape 

of  the 

£.e. 

The  singularity 

of 

Order  of 

for 

F  -*•  0 

Order  of 

a2  foi 

F  ■*  0 

Order  of 

ratio 

Order  of  the 

straining 

u^£ 

(z-l-ih)"1 

e-^2 

c/F2 

r. 

/z  -  1  -  ih 

e_h/F2 

t/F 

eF 

> 

In  (z-l-ih) 

mm 

e~h/FVln2F 

e  In  F 

G  T  In  F 

jkfri!"*  . . . 
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The  column  before  the  last  summarizes  the  main  findings:  the 
quantity  appearing  there  has  to  be  small  in  order  to  ensure  that 
the  usual  linearized  thin  body  approximation  is  uniform.  It  is 
worthwhile  to  mention  that  in  all  the  examples  in  which  detailed 
computations  have  been  carried  out  so  far  (Tuck,  1965,  for  a 
circular  cylinder;  salvesen,  1969,  for  a  hydrofoil  and  Dagan, 
1972a  for  a  source),  the  shapes  were  blunt. 


5 .  Derivation  of  Uniform  Small  Froude  Number  Solutions 

(i)  Valoait ^  etrainin g 

The  nonuniformity  of  the  free  waves  expansion  is  assumed  to 
originate  from  the  illegitimate  expansion  in  an  r.  power 
series  of  an  exponential  of  type  exp  [ icF  ( z)  /F2 )  for  e/F2  ■>  o  ( 1  )  : 
Such  a  term  may  result  from  the  straining  of  the  free-surface 
velocity  due  to  the  presence  of  the  body.  To  make  the  idea  more 
precise  let  us  replace  the  first  order  free  surface  condition  (1.2) 
by 


Im  [iF2  ( 1-t-tA)  -  f]  *  0  (y  =  0  )  (1.33) 

where  1 + eA  is  a  strained  velocity.  A(z;h,F)  is  assumed  to  be 
analytical  and  A 0  for  x  -*• «  .  A  is  obviously  related  to 
and  under  the  usual  thin  body  expansion  (1.1)  it  is  cast  in 

(x)  (1.6).  vv'e  now  keep  it  in  the  first  order  equation  and  ; 

solve  for  f  with  boundary  conditions  (1.4),  (1.5)  (again,  we  ; 

! 

do  not  consider  the  second  order  body  effect  (1.9)).  j 

i 


t 

| 
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The  solution  is  immediately  obtained  as  follows 


„  rj.Z  .  if  B(o)  ,  i  [  du 

£  -  +  /  Y+^XToT  e  p(~  F1*  !  r+TXm 


du 


)  do  3 


(1. 39) 


where 


B  (o) 


£  £. 

*  (w.+w^)  (1+r.Re  A)  ">  (w.+w^)He  A 

/  — L~i— - —  dp  =  2 irt  w“(z)  +  e  /  — -  dP  (1.40) 

a.  oo  U  M  J.  _  OO 


o  -  p 


For  e  =  o(i)  we  are  generally  entitled  to  expand  in  a  power 
series  the  terms  B(A)  and  (1  + eA)  ^  in  (1.39)  and  (1.40). 

This  yields 


■f  - 


If  Ji  fj  J 


i  X  /T2 


exp(ff  } 

o 


ACu)dy)daJ  +0(t!)  (1.41) 


We  did  not.  expand,  however,  the  lart  exponential  in  (1,39)  since 
under  the  limit  F  -*■  0  the  ratio  e/F2  is  not  necessarily  small. 
£q.  (1.41)  proves  our  assertion  on  the  effect  of  a  first  order 

velocity  straining. 

The  uniformization  procedures  of  the  small  Froude  number 

solution  suggested  in  previous  works  are  underlain  by  similar 

ideas.  Ogilvie  (1968)  has  arrived  to  a  free-surfece  condition 

similar  to  (1.38)  by  intuitive  reasoning:  as  the  wave  length  of 

th v.  free  waves  becomes  small  compared  to  the  length  scale  of  the 

velocity  field  ,  the  velocity  variation  has  to  be  included 

in  a  first  order  approximation.  Moreover,  it  was  suggested  that 

£  u 

os  F  0  A  =  2  (w.  +  w.)  for  y  «  0 ,  since  for  F  =  0  w.  degene- 

l  u  A 

rates  into  the  rigid  wall  solution  wicwi+wi‘  An  equation 
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6itnilar  to  (1.38)  is  then  obtained  from  the  exact  free-ourface 
condition  by  expanding  with  l+c(w^+w^)  as  the  basic  unperturbed 
velocity  field. 

Dagan  (1972b)  has  arrived  at  a  similar  result  by  using  a 
quasi-linear  equation  as  a  model  of  the  nonlinear  f ree-surf ace 
condition  (in  both  works  the  more  ambitious  task  of  solving  the 
problem  of  a  small  Froude  number  flow  past  a  body  of  finite 
thickness  has  been  undertaken) . 

Although  the  arguments  are  plausible  in  principle,  the 
assumption  that  A  =  z(w^-»-w^)  is  opan  to  criticism,  since  it 
has  been  shown  here  that  the  degeneracy  into  the  rigid  wall 
solution  is  not  uniform. 

Instead,  we  are  go-ng  to  determine  here  the  straining 
function  A(z)  in  a  different  way.  We  assume  that  for  x  ■+■  -  “ 
f  in  (1.41)  includes  the  first  order  term  f^w  (1.36),  as  well 
as  the  lowest  order  term  in  F  appearing  in  £_,w  (1.23).  We 

require,  therefore,  that  under  an  additional  expansion  for 
eA/F*  <=o(l),  (1.41)  should  degenerate  into  the  thin  body  expan¬ 

sion  (1.1).  Expanding  the  exponential  in  (1.41)  gives  for  x  *  ->*> 


,.w 


:e  e 


-iz/F‘ 


/  w^  (o)e^0/'F  da  - 


2ie2  -iz/F2 

T1-* 


oo 

/  A(o)dc 

—  to 


/ 


it/F2 


dt  +  0(e2) 


(1.42) 


The  first  term  in  (1.41)  recovers  the  first  order  solution  e f ^ 

(1.36).  The  second  term  may  be  identified  with  the  F  lowest 
w 

order  term  of  f^  (1.25),  which  may  be  written  as 


-a 


J 


I 

4 


i 


! 

1 t 
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e’f*  =  c2e'i2/F2  f  Uw^(c)]2eia/F2  - 


-  pr  fw“(o)  +wj(cj]  J  w^(p)eip/F  dp}do 


ip/F- 


(1.43) 


Identification  of  A  such  that  the  integrands, 
F  term  of  P^Co),  in  (1.42)  and  (1.43),  become 


.f2  (w?)2  eiz/F  t  u 

A  ( z ;  h  ,  F )  =  ±4 - - -  +  2(w.V)  = 

^  i  u  ix /F2  .  1  1 

J  w.  e  dr 


(w”)2 

~u  I 
w.  w,  +  w, 

J.  J  1 


+  2  (w*+w^) 


i.e . ,  the  lowest 
identical  yields 


(1.44) 


Hence,  the  velocity  straining  function  A  i6  found  to  be 

diff  rent  than  that  suggested  in  previous  works,  which  included 

only  the  last  term  of  (1.44),  but  failed  to  take  into  account 

9  u 

the  singular  term  related  to  (w^-w^)  in  (1.25).  The  reasoi 
is  quite  transparent:  this  last  term  is  identically  zero  for 
y  =  0  according  to  the  rigid  wall  sol.ution.  The  rigid  wall 

solution  is  not  uniform,  however,  and  the  singularities  of  f. 

£  u  ^ 

(1.17)  and  of  f ^ 4  f^  are  different  at  the  location  of  the 

image  of  the  body  across  the  f ree-surf ace . 

Obviously,  for  F  -*•  0  and  cA/F2  *»  0(1)  the  exponential 
in  (1.41)  cannot  be  expanded  like  in  (1,42).  By  the  same  token, 
the  transfer  of  the  velocity  straining  factor  eA  into  the 
right  hand  side  of  the  second  order  free -surface  condition  (1.6) 
is  not  legitimate  for  small  Iroude  numbers  'n  general. 
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Computing  the  wave  amplitude  with  the  aid  of  (1,41)  is  a 
difficult  task  in  the  two-dimensional  case  and  becomes  extremely 
tedious  in  three-dimensions.  For  these  reason  we  consider  sub¬ 
sequently  a  simplified  procedure  for  rendering  the  solution  uniform. 


(ii)  Coordinate  straining 

We  assume  now  that  the  exponential  terms  which  cause  the 
small  Froude  number  nonuni formity  are  a  result  of  a  coordinate 
straining.  Lighthill  method  (3ee  Van  Dyke,  1964)  implies 
an  infinitessimal  straining  of  the  physical  plane  and  deriving 
the  straining  function  from  the  equations  of  flow.  We  adopt 
here  a  modified  technique  applied  by  Van  Dyke  (1964)  to  the  case 
of  inviscid  flow  past  airfoils:  we  carry  out  the  straining  in 
the  solution,  rather  than  in  the  equations,  and  determine  the 
straining  function  from  the  requirement  that  the  second  order 
term  should  not  be  more  singular  than  the  first. 


To  illustrate  the  method  we  begin  with  the  example  of  a 
source-sink  body  (Section  1,3).  The  straining  has  as  effect  a 
virtual  displacement  of  the  images  of  the  two  singularities  from 
?.£,  z^  to  +  6  ,  zt  +  6  zt  ,  respectively,  with  6z^  =  0(e) 

and  <5zt  =  0(e), 

The  first  order  term  of  the  free  wave  expansion  becomes  now 
by  using  (1.17)  and  (1,27)  for  x  -*• 


f 


w 


->  00 

e~lz/F*  j  , - 1 - 


o-z 


io/F2 

e 


do 


(1.45) 


For  e  =  o(l)  and  F  fixed  we  can  expand  in  (1.45)  and  obtain 
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e  -iz/F2  /  ,  1  1  .  io/F2  , 

-  e  )  (* — - - — — )  e  do  + 

a-Z|  cr-zt 

v  -iz/F  /“,«*£  _  ,  ic/F: 


/  - r^— )  e 

-®  (o-z^)2  {o-zt)J 


fw  2e  -iz/F2  -  1Z£/F2  --  1Zt/F\ 

^  ■  p  e  (&zi  e  -  Szt  e  ) 


(1.46) 


Hence,  the  first  term  of  (1.46)  recovers  f^w  (1.31). 
Consequently,  the  second  order  term  of  the  free  waves  potential 
will  be  made  up  this  time  from  f™  (1.32)  plus  the  last  term 
of  (1.46),  provided  that  the  straining  is  of  order  £  . 

We  determine  now  and  <5zt  from  the  requirement  that 

the  term  of  order  e/F2  in  the  amplitude  of  the  free  waves, 
which  is  the  origin  of  the  small  F  nonuniformity,  should  vanish 
in  the  solution,  separately  for  the  source  and  the  sink.  We  thus 
obtain 


27  {b2  +ic2} 


(1.47) 


27  (b2  "  ic2) 


(1.43) 


where  t>2 ,  c2  are  given  in  (1.32). 

The  uniform  first  order  solution,  valid  for  e/F2  =  0(1), 
is  easily  derived  from  (1,45) 


w  -iz/F2  1  ) /F 

r  =  2i  e  1Z/t  [e  *  * 


i  (2  + 6z  )/F* 

-  e  *  c  ]  (1.49) 


By  using  (1.47)  and  (1.48)  we  finally  obtain  from  (1.49) 


-  [h-  (c9e/2tt)  ]/f2 
fw  =  -1e  * 


l-(cb2/2r)  2 

sm  - rn -  e 


(1.50) 


.  ,  uumrnmmmMnitPtyW?' (S'11'1' 
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The  coefficient  c2  (1.32)  is  associated  with  the  second 
order  interaction  between  the  source  and  the  sink;  it  has  the 
effect  of  diminishing  the  effective  submergence  depth  of  the 
body  {when  h+0  this  term  vanishes).  bn  is  associated  with 
the  nonlinear  effects  of  the  leading  and  trailing  edges  singulari¬ 
ties  upon  themselves.  It  manifests  in  an  apparent  change  cf  the 
body  length  and  consequently  in  a  shift  of  the  curve  of  the  ampli¬ 
tude  (and  wave  resistance)  as  function  of  F. 

Again,  (1.50)  shows  that  the  small  Froude  number  non uniformity 
is  a  result  of  an  illegitimate  expansion  of  the  exponential  and 
trigonometric  functions  in  (1.50)  in  an  e  power  series  for 
e/F2  «=  0(1).  In  other  words,  when  the  straining  becomes  of  the 
order  of  the  wave  length,  it  has  to  be  maintained  in  the  first 
order  approximation. 

we  are  going  now  tc  generalize  the  procedure  for  a  body  of 
arbitrary  thickness  distribution  t  (x  )  «=  dt/dx.  The  straining 
has  now  as  effect  a  continuous  inf initessimal  displacement  of 
the  image  of  the  body  skeleton  from  z  to  z  + 5z  .  The  first 

8  8  S 

order  solution  ((1.16)  and  (1.17))  becomes 


*w  2e  -iz/F2  . 

f  *=  -  — -  e  I  do  e 

TT  * 


ia/FJ 


/ 

-1 


dx_ 


T(V 


a  -  2  -  6 z. 


(1.51) 


where  z  =  x  +ih.  The  straining  has  been  taken  into  account 
only  in  the  denominator  of  (1.51),  because  only  the  residues  at 
c  -  z  -  6z  ^  0  are  contributing  to  the  lowest  order  terms  in  F. 
An  inf initossiroal  change  of  the  limits  of  integration  or  of 
t (x  )  in  (1.15)  yields  higher  order  terms. 
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For  6t  «*  0(e)  and  F  fixed  the  integrand  in  (1.51)  may 

8 

be  expanded  as  follows 


-W 

f  • 


if  e-^rl  l" do  /  dx  -i—  - 

n  Lis- 

•■00  —  ^ 


a  -  x. 


2e  -iz/F 

T*  e 


2  00  1  t6  z 


f  da  f  d  x 


-1 


S  (a  -  i  )2 

B 


(1.52) 


-4ic  e-i2/lJ  / 1  T(X  le^d*  +4|e-i2^!  /  ,(* ,  >«5 
^2  s  s  r  ^2  s 


-  izs/F 
e  s  dx 

8  b 


.  j 


The  first  term  in  (1.52)  is  precisely  f^w  (1.36);  the  unknown 
straining  function  62g  is  now  determined  from  the  requirement 
that  the  second  term  of  (1.52)  should  cancel  the  lowest  F  term 
of  f .,  (1.43) ,  i.e,( 

1  _  i*a/F* 

_{  T'Xs'SZ3  6  dX8  “ 


(wj+w“)  j  w“  u(~f>  dp)e10/F2do  (1.53) 
•■00 


i 


I 


To  determine  <$z  in  a  simple  way,  advantage  is  taken  of  the 
fact  that  the  lowest  F  terms  in  (1.53)  are  associated  with  the 
singularities  of  the  edges  (an  intermediate  point  of  discontinuity 
can  be  easily  accounted).  What  matters,  therefore,  is  6z^  and 
6z^  .  Any  continuous  straining  between  the  edges  is  acceptable 
as  far  as  the  most  singular  terms  are  concerned.  Assuming,  for 
the  sake  of  simplicity,  a  linear  straining  we  have 


s\-izt  .  ezi*s\ 

- 3 - xs  +  - 3 - 


(1.54) 
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Substitution  of  (1.54)  into  the  left  hand  side  of  (1.53) 
yields  for  the  lowest  F  term 


1  i2g/F2 

(5z.  +5z  )  /  t  ( x  )  e  dx  . 

v  ^  S  D 


(1.55) 


Equating  the  lowest  F  terms  resulting  from  the  integration  in 
the  r.h.s.  of  (1.54)  and  from  (1.55)  renders  in  an  unique  manner 
+  .  An  additional  relationship  is  obtained  from  the 

requirement  of  separate  cancelation  of  the  leading  and  trailing 
edges  waves  (obviously,  for  a  fine  shape  of  the  trailing  edge 
6zfc  *=  0)  . 

The  estimates  of  Section  1.4  permit  to  evaluate  the  order 
of  the  straining  6z^  for  different  types  of  leading  edge 
singularities.  The  results  are  given  in  the  last  column  of 
Table  1.  The  straining  becomes  weak  as  the  shape  of  the  leading 
edge  becomes  fine  and  it  is  F  dependent,  excepting  the  source¬ 
like  case. 


6 .  Conclusions 

It  has  been  shown  that  the  slenderness  small  parameter  e 
appears  in  the  expression  of  the  potential  of  the  free  waves 
generated  by  a  submerged  body  not  only  in  the  amplitude,  but 
also  as  the  ratio  e/F2  in  the  wave  number.  Like  in  other 
problems  characterized  by  two  scales  (Cole,  1968)  a  power 
expansion  in  e  does  not  render  a  uniform  solution  unless 
e/F2  **  o(l);  this  last  estimate  has  been  sharpened  and  shown 
to  depend  the  nature  of  the  leading  (and  trailing)  edge 
singularity. 
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Two  procedures  of  rendering  the  email  Froude  number  solution 
uniform  have  been  suggested;  f ree-surf ace  velocity  straining  and 
coordinate  straining.  The  first  procedure  has  the  advantage  of 
making  uniform  the  second  order  preesure  term,  whose  integration 
yields  the  amplitude  of  the  free  waves;  for  this  reason  the 
velocity  straining  is  easily  expressed  with  the  aid  of  the  first 
order  solution.  Computing  the  uniform  solution  is,  however, 
extremely  difficult.  The  coordinate  straining  ensures  the  uni¬ 
formity  of  the  expansion  of  the  free  wave  amplitude  (and  the 
wave  resistance)  and  the  computation  of  the  straining  factors  is 
more  difficult  than  in  the  first  case.  Once  determined,  however, 
they  provide  immediately  the  uniform  first  order  solution. 

Two  problems  have  not  been  touched:  (i)  the  second  order 
body  effect, and  (ii)  circulation.  As  for  (i)  it  has  been  shown 
previously  (Dugan,  1372a)  that  the  body  correction  is  uniform  as 
F  -*•  0  (e.g.  ,  for  a  source  e2a2  c  °^e2  )•  Circulation  may 

have  an  important  influence  on  the  wave  amplitude  (Salvesen,  1969). 
It  can  be  treated  similarly  to  the  thickness  with  no  problems  of 
principle.  We  have  purposely  considered  the  effect  of  thickness 
solely  because  the  two-dimensional  solution  serves  here  only  as 
a  case  study  for  the  three-dimensional  flow  problem.  In  appli¬ 
cations,  two-dimensional  flows  are  generally  at  high  Froude  numbers. 
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THREE-DIMENSIONAL  FLOWS 


1.  The  Thin  Body  Expansion  and  Fourier  Transforms 

Let  now  z  be  a  vertical  coordinate ,  while  the  axis  x 
and  y  lie  in  the  horizontal  plan**  of  the  unperturbed  free 
surface.  ia  the  dimensionless  velocity'  potential  which  is 

expanded  in  a  thin  body  approximation  as  follows 

4>(x,y,z;e,F ,...)  =  -x  +  e^fx/y/^F,,..)  +  £2<f)2(x,y,z;F,...)  +  ...  (2.1) 

where  F  *  u'/(gL')  ,c  is  again  the  length  Froude  number. 

The  free  surface  conditions  satisfied  by  the  harmonic 
functions  and  <t>0,  given  here  for  the  sake  of  completeness, 


3  <t>^  3  $  ^  3  $  ^ 

F  +  JT  ■  y  ” 


(z«0)  (2.2) 


32$2  ?.<l>2 

+  TT  -  *  3x~  3  P2(x^} 


dU.  dV  dV  dw.  dMi. 

-  -F*(3ul  5x  +  2vl  TT-+U i  7r+2“i  51ir'p!ul 

where  u^,  vx#  wi  are  t*ie  ^^rst  or<3er  x,y,z  velocity  components  and 
u  0  represents  the  "artificial  viscosity"  added  in  order  to 


. . . . . 


HYDRONAUTICS ,  Incorporated 


25. 


satisfy  the  radiation  condition  (we  introduce  it  rather  for  an  easy 
account  of  the  integration  paths  in  the  complex  plane) . 

The  body  boundary  conditions  differ  depending  on  whether  one 
considers  slender  or  thin,  submerged  or  surface  piercing  bodies. 

For  this  reason  they  will  be  formulated  separately  in  the  sub¬ 
sequent  examples. 

Again,  like  in  the  two-dimensional  case,  we  consider  here 
the  part  of  the  second  order  solution  which  is  regular  in  z  <  0 
and  satisfies  (2.3),  and  do  not  investigate  the  second  order  body 
correction  as  well  as  the  line  integral  correction,  since  the 
nonlinearity  of  the  free  surface  condition  is  the  most  severe 
as  F  0  . 

We  summarize  now  the  different  Fourier  transforms  to  be  used 
in  the  sequel;  with  ®  denoting  the  Fourier  transform  of  <f> 
with  respect  to  x,y,  we  have 


<j)  (a,  B ,  z)  »  ~  /  dx  /  dy  <p  (x,y ,  z)  e1  (ox+ey }  /F 


(2.4) 


<Mx,y,z)  =  -p-  /  da  /  dB$  (a, 8 , z ) e"* (ax+By } /F2  (2.5) 

mOO  —00 

where  a, 6  are  coordinates  in  the  transform  plane.  We  shall  also 
use  the  convolution  transform  which  may  be  written  as 

, _ _ _ ^  i  00  “ 

<}>  (x,y)  (x,y)  -  /  dv  /  dx  4- (v ,  t)  ^  (u-v ,  B-x)  (2.6) 


where  <p,ip  are  arbitrary  functions  which  have  FT  (Fourier 
transform) . 
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Let  now  ♦  *  represent  the  first  order  potential  of  flow 
past  the  body  in  an  infinite  domain  and  the  same  for  the 

flow  past  the  image  of  the  body  across  the  plane  z«0.  Then, 
the  well  known  (Wehausen  and  Laitone,  1960)  first  orde-  solution 
$  .  satisfying  Laplace  equation  and  (2.2)  becomes 

♦l  *  +  ♦!  (2.7) 


where 


=  -2  F2  — 

1  -2. 


;?(«.#. o)«pi/p! 


(z<0)  ,  (2.8) 


a*-  p  -  ipa 

W1  B  "nd  P  =  6 ' .  By  (2.5)  we  have 

,r  1  ,  f°° w1(a,6,0)  -i  (ux+gy+ipz)/F2 

y,  -  -  — yrr  »  uu  i  ap  —  e 

i  it*  _i  a  P“lva 


(2.9) 


being  understood  that  u  -*■  0  in  the  final  expressions  ($£  is 
regular  for  z<0  and  0  for  z  ■*■*•”)  . 

The  FT  of  the  second  order  solution,  regular  for  z  <0  and 
satisfying  (2.3),  is  similarly  given  by 

FJP2(a,B)  epz/r' 

02(tt,s.z) - yr -.-p  i-iua -  <2-10! 


From  (2.3)  and  (2.6)  P^(o,6)  may  be  written  after  a  few 

integrations  by  parts  aB  follows 

oo  00 

-  YttF*  I  /  dx  {  [  3  (a-v  )  v2  +  2  (6-x ) vx  +  (a-v)  x2  - 

•  co  -»oo 

-  2 (a-v) 2  v  3  -  (a-v) v4] (v,x ,0) (a-v,B-x,0)}  (2.11) 


P2 (a, 8) 
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It  seems  that  all  the  terms  of  P2(x,y)  (2.3)  contribute  at  the 

same  F  order  in  the  expression  of  Pn(a,B).  It  is  sufficient, 
therefore,  to  single  out  one  of  these  terms  in  order  to  estab¬ 
lish  the  asymptotic  behavior  of  $  as  F^O.  Like  in  the  case 
of  two-dimensional  flow  (Sec.  15),  it  is  convenient  to  select  the 
term  w,  (3w^/3x)  and  to  write 


3w, 


P2(x,y,F)  ^  -2Fltw1 


(z  =  0)  (2.12) 


being  understood  that  (2.12)  expresses  the  asymptotic  dependence 
of  P0  on  the  variable  F. 

Finally,  it  is  sometimes  convenient  to  operate  with  polar 
coordinates  in  the  transform  plane,  a  »  p  cos  9  and  B  •=  p  sin  0  . 
Then,  if  the  potential  <fiW  =  +  eH”  of  the  far  free  waves  is 

written  in  the  following  form 

-it/2 

=  i  Rs  J  de  [Q(e)+iF(one-1(xsecety  sln«»ecJe+iz  (X 

-,/2  (2-13> 

the  dimensionless  wave  resistance  is  given  by  (Maruo,  1966) 

-tt/2 

R-  - — -  *  I*  I  [Q2  (6)  +  P2  (6)  ]secs9  d9  (2.14) 

0 . 5p 1 U*  *L' 2  n  Jn/ 2 


where  p'  is  the  fluid  density.  The  amplitude  functions  P  and  Q 
result  froi 
written  as 


result  from  the  thin  body  approximations  and  $2  and  can  be 


P  «  eP.  +  c2P„ 
1  ^ 


Q  =  +  e-Q2 


(2.15) 


I . . 
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which  leads  by  substitution  in  (2.14)  to 

R  =  +  e*R2  +  . . .  .  (2.16) 


2.  Small  Froude  Number  Solution  for  An  Isolated  Source 


(i)  First  order  solution 

Like  in  the  two-dimensional  case  we  begin  the  study  of  the 
limit  F  -*•  0  of  the  thin  body  solution  with  the  example  of  an 
isolated  source,  because  of  its  simplicity  and  because  the  source 
is  the  fundamental  singularity  underlying  any  slender  or  thin 
body  solution. 

The  potentials  of  flow  in  infinite  domains  are 


4> 


l 

Is 


1  _ 1_ _ 

[x2  +  y2  +  (z+h) 2 


(2.17) 


u  _  _  i 
Is  Tn 


_ _ 1 _ 

(X2  +  yJ  +  (z-h) 2  ] 


(2.18) 


The  source  is  located  at  x  =  y  =  0,  z  = -h  (the  reference  length 
L'  is  left  unspecified) .  In  the  frame  of  the  first  order  approxi¬ 
mation  <p  =  -x  +  represents  the  flow  paet  a  slender  body  of 

revolution  with  a  blunt  nose  and  of  semi-infinite  length 
(-®°<x<0).  The  small  parameter  c  is  equal  to  TT(r'/L')2,  where 
r'  is  the  radius  of  the  body  circular  cross-section. 


We  consider  now  the  complex  a  plane.  The  function 

p  =  /a*  +  has  its  branch  lines  depicted  in  Fig.  2  with  the 

values  of  the  argument  selected  such  that  for  a  ,  6  reals  p  is 

real  and  positive,  such  that  the  condition  ■*  0  for  z  ■*  -  « 

is 

is  identically  satisfied. 
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The  poles  ir  (2.20)  are  a  =*  tn  3  y 
end  have  to  be  circumvented  from  below.  We  can  now  awing  the 
integration  path  from  the  real  axis  in  tho  a  plane  to  the 
branch  cuts.  The  result  is  easily  found  to  be  made  up  of  the 
so  called  local  disturbance  loc  and  the  free  waves  term 

♦is 


♦  /  e"B^|/F*  cos(8y/F2)de 


f°°  e  9  { (S+s)  2Bin[y(z-h) /F2 )  +  v  oos  (y  (z-h)/F2 ] } 


(6+8)“  +  y2 


(2.21) 


where  y  =  /s* +  2  8s  ,  and 


a2 (z-h)/F 

„  «=>  sin  (a  x/F2)co6(6y/F2H5  p 

f  — - _ _ 

w  is  ttE“  '  =T 


2a  -  a 
F  F 


=  |  sin [  (x  oos  6  ty  sin  0  sec2e)/F2]e(z“h)sec26/^2sec26  de 

-ir/2 

(x  <  0)  (2.22) 


Like  in  the  two  dimensional  case  ^  ^oc  is  regular  for  z  <h 


and  tends  to  zero  algebraically  for 


After  a  change  of 


variables  it  can  be  expanded  uniformly  in  an  F2  power  series  for 

F  0  and  fixed  x,y,z  (with  x2  +  y 2  +  ( z-h)  2  >  0);  as  a  result 
,r  loc  _.u  _  . 

♦is  "  2*l<»  +  0(F  >* 

w 

The  free  wave  potential  <f>^3  is  different  of  zero  only  for 
x<  0.  BecauBe  of  the  exponential  term,  it  is  negligible  with 
respect  to  4> oc  excepting  in  an  arbitrarily  thin  wedge  z-h  --  6x 
(x  <  0,  6  arbitrarily  small) .  Again,  like  in  the  two-dimensional 
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case  $*g°C  i*5  not  ah  uniform  approximation  of  <f>^g  and  for  x  -*■  -< 

dominates  the  solution. 


£  li 

Summarizing,  as  F  -*  0  <J>^3  =  +  0(F2),  where 

+  is  the  rigid  wall  solution,  but  the  limit  is  not  uniform 

for  z-h  =  fix  (x  <  C)  ,  where  it  has  to  be  supplemented  by  <{>Y  . 


(ii)  Second  order  eolation  (the  free  v>~.v*e  potential) 

In  contrast  with  the  two-dimensional  case,  now  it  is  not 
possible  to  obtain  the  solution  in  a  closed  form.  It  is  relatively 
easy,  however,  to  evaluate  the  order  of  <f>2  for  F  -  0  by  retaining 
only  part  of  the  terms  of  P2  (2.12),  as  it  has  been  done  in 
Section  14.  Like  ir»  two-dimensions  w^  may  be  written  as  (see  (2.7)) 


w,  -  w. 
Is  Is 


(2.23) 


where  in  the  case  of  a  source  (2.17-2.18) 


z  ♦  h 


(x2  +  y2  +  (z+h) 2 )  3/' 


(2.24) 


z  -  h 


(x2  +  y 2  +  (z-h)1] 


(2.25) 


and  wf_  =  34>5\/9z  (2.20).  After  substituting  w.  in  P„(x,y) 

(2.22)  it  can  be  shown  that  the  term  -2F2w!r  (3w,  /3x)  contributes 

IS  J.S 

to  the  lowest  order  term  of  <P2  as  F  0 .  Hence,  by  (2,22)  and 
(2.25), 
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3w 


u 


?2(a,6)  -  -2F^b  -jit  -  4ic(w“s)> 


(2  =  0)  (2.26) 


$2(a,8)  ^ 


-4ia  F2  [w^s  (a.,  8)  ) 2  4iF2  [w^g  (p  ,  6 )  ) 2 sec  6 


a  - p - iya 


p -  sec7  6  -  iu  sec  0 


(2.27) 


The  Fourier  transform  of  (w^g)2  ea£ily  found  in  polar 

p,0  coordinates.  From  (2.2  5)  we  have 


u 


(w'{B<x,y,0>)!  = 


2  \  3 


(x2  +  y2  +  h2 ) 


and  by  (2.4) 


,  .2  tt  oo  ior  cos(X-9)/F2 

<“u>!  -  357T  /  a»  I  r° 


dr 


(r2  t  h2 )  2 


(2.28) 


(2.29) 


where  x  =  r  cos  X  ,  y  =  r  sin  X  .  Integration  in  (2.29) 
(Gradshteyn  &  Pyzhik,  1965)  yields 


,  .u  x2  _  h*  ;  ,  ,or,  r  dr  p2  ,phx  /0 

(wis)  "  TTP  J  Jo  f^*  '  HW  K-2  iN*  (2*30) 

o 


where  J  and  K  are  Bessel  functions  of  the  first  kind  and 
o  -2 

of  the  second  kind  (modified),  respectively.  By  using  (2.27), 
(2.30)  and  (2.5)  we  obtain 

.  .*  -  p>K_2(ph/ti)e‘ip<X0“®*»8in9+i*)/f2Mc9 

^2s  u  6  4tj  5F6  /  ^  dp  ”  ~  sec7$" -  Tu  secTS 


Like  in  the  two-dimensional  case  it  can  be  shown  that  the  near 
field  expansion  of  (2.31)  yields  a  uniform  F2  power  series  expan¬ 
sion  of  <fc 

8 


.31) 
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The 


free  waves  term  of  (2.31)  may  be  written  ae 

4  & 


tr/2 


4>2g  x  YgffTgr  /  sec®0K_2(h  sec20/F2  )coe  [  (xsec0+ysec0sec20)/F2  ] 
-tt/2 


x  ez  sec20/F23ecQde  (2.32) 


As  F -*•(),  K_2  (h  sec28/F2)  can  be  expanded  in  an  asymptotic  series 

which  yields  for  the  leading  term  of  <p^s 

♦  oe  %  ^  /  sec60  e cx)6 [ (xsec0+ysin0sec20)/F2]d6  +0(=rr)  .  (2.33) 
zs  *  -n/2  * 


(Hi)  UitoueBion^o^^reBultB 

We  are  now  in  a  position  to  discuss  the  final,  and  most 
important,  topic  of  this  section,  namely  that  of  the  uniformity 
of  the  free  waves  expansion  <f>w  =  -+•  c 2 4>^  under  the  small 

Froude  number  limit.  Eqs.  (2.22)  and  (2.33)  show  that  the  thin 
body  solution  is  not  uniform  unless  e/F3  =  o(l);  the  latter 
condition  originates  from  the  comparison  of  the  amplitude 
functions  P^(0)  and  Q2(8)  in  (2.22)  and  (2,33),  respectively. 
The  same  condition  ensures  the  convergence  of  the  coefficient  of 
wave  resistance  (2,16).  Taking  the  submergence  depth  as  reference 
length,  we  can  write  in  term  of  variables  with  dimensions 


U: 


r ’  2g2 

7^7171 


(2.34) 


where  r’  is  the  radius  of  the  body  generated  by  the  source. 
This  ratio  has  to  be  compared  with  the  criterion  T'g/U’2  »  o(l) 
which  ensures  the  uniformity  of  the  thin  body  solution  in  two 
dimensions. 
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We  can  continue  and  analyze,  like  in  Section  14,  uniformity 
criteria  of  the  thin  body  expansion  for  F  -*•  0  in  the  case  of 
bodies  of  revolution  with  milder  singularities  at  the  nose  (and 
the  tail) .  From  the  applications  point  of  view,  however,  it  is 
of  interest  to  focus  the  analysis  on  the  case  of  thin  ships  rather 
them  on  '..hat  of  submerged  slender  bodies. 


3 .  Small  Froude  Number  Solution  for  Thin  Ships 
(i)  General.  First  order  solution 

Let  y=  ±cf(x,z)  be  the  equation  of  the  surface  of  the  ship, 
where  c  =  B'/L'  is  the  beam  length  ratio.  The  first  order 
velocity  potential  <J>1  is  obtained  by  integration,  over  the 
area  of  the  center  plane  S  ,  of  the  potential  of  an  isolated 
source  <t>ls 


a  =  -2  /  /  <fr  (x-x,y,z-z)  dxdz  (2.35) 

1  S  dX  19 

where  4>ls  =  4>^s  -  4>^g  +  is  given  in  (2.17),  (2.18)  and  (2.19) 

with  z  replacing  -h.  For  the  sake  of  simplicity  we  consider  a 
symmetrical  ship  solely. 


The  essential  difference,  in  the  present  context,  between  <J>. 


in  all  the  other  cases  (of  submerged  bodies) 


(2.35)  and 

considered  in  the  preceding  sections,  stems  from  the  fact  that 
the  body  is  now  piercing  the  f ree-surface .  Consequently,  the 
exponential  factor  e  and  e  ^  sec  ,  present  in  the 

two-  and  three-dimensional  solutions,  respectively,  does  not 
appear  anymore  in  (2.35)  for  the  waterline  (z  «=  0)  singularities. 
For  F 0  and  h  fixed,  the  exponential  factor  ensured  previously 
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the  rapid  decay  of  the  free  waves  amplitude  and  of  the  wave 
resistance  coefficient  no  matter  how  blunt  the  shape  was.  It 
also  ensured  the  separation  of  the  near  field  expansion  and  of 
the  free  waves  potential. 

In  the  case  of  a  thin  ship  the  shape  has  to  be  sufficiently 
fine  in  order  to  ensure  the  finiteness  of  the  velocity  potential. 
It  can  be  shown,  for  instance,  that  for  a  source-like  bow  shape 
the  second  order  potential  is  not  integrable.  Since  usual  shapes 
are  far  from  being  so  blunt,  we  shall  limit  the  discussion  here 
to  wedge  like  bodies,  i.e.,  with  singularities  associated  with 
finite  entrance  and  shoulder  angles  at  most  (the  case  of  an 
elliptical  bow  is  also  less  interesting) . 

Further  simplifications  of  the  analysis  arc  achieved  if  we 
take  into  consideration  the  well  known  asymptotic  properties  of 
4^  (2.35)  as  F  -*•  0  (see,  for  instance,  Lunde ,  1963): 

(i)  the  dominant  contribution  originates  from  the  shape  at  the 
waterline  f(x,0),  and  (ii)  from  the  singularities  of  f(x,0)  , 
i.e.,  from  the  points  of  discontinuity  of  3f(x,0)/3x.  For  this 
reason  we  consider  the  simple  example  of  flow  past  a  wedge-shape 
cylindrical  bow,  i.e., 

x  (-1  <  x  <  0 ,  -h  <  z  <  0) 

(2.36) 

(  X  <  -1,  -h  <  2  <  0) 

the  reference  length  being  the  forebody  length.  f  (2.36) 
incorporates  the  essential  features  of  any  smooth  shape  between 
x  ■  0  and  x  =  -l  with  the  same  angle  discontinuities.  The 
influence  of  the  stem  is  also  disregarded  because  as  will  be 
shown  later,  the  most  singular  terms  are  associated  with  the 


f(x,z)  = 
f(x,z)  = 
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interaction  between  the  singularity  at  x»0  (or  x  » -1)  with 
itself  and  not  with  the  interaction  between  the  different  singu¬ 
larities. 


From  (2.35),  (2.19)  and  (2.36)  we  immediately  obtain  for 
♦j  -  *i  - 

>1  -  - 

■  -  £  / «  /“ j^^^ccose/F2' 


-IT 


x  e“ip  (x  006  6  +y  sin  6  +iz)/T2secJe 


(2.37) 


The  near  field  expansion,  i.e.,  the  local  disturbance  A*  *oc  may 
be  found  like  in  (2.21)  and  again  V^'00  =  2#^  +  CHF2).  The  free 
waves  potential  is  found  from  (2.37)  as  follows 


♦ 


w 

1 


si*  r 

n  -rr/2 


(l-e'h  sec‘^)  (i-e1  ^/F2)  , 

v  i  (x  oos  6  +y  sin  0  -iz)/F2 
x  e  J  ooe 


(2.38) 


As  F  0  a  lowest  order  term  in  (2.38)  is  obtained,  for 

instance,  from  integration  over  e^x  cos  ®  s^n  ®  COB  q  f 

i.e.,  from  the  potential  of  a  wedge  of  infinite  draft  and  of  con- 

stant  aperture  angle,  which  leads  to  <=  0(F*)  for  fixed 

y ,y  (x2+y2  /  0)  and  z=0.  By  (2.14)  the  coefficient  of  wave 

resistance  derived  from  this  term  is,  at  first  order,  R  =  0(e2F4). 

w 

Hence,  for  z  »  0  <f> ^  is  no  more  exponentially  small,  as  F—  0 , 

in  comparison  with  the  near  field  solution;  moreover,  for 
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(x2+y2)/F2  ■  0(1)  and  z*0  it  becomes  of  the  same  order  as  the 
second  terra  of  the  near  field  solution,  i,e.  ,  0(F2).  For  any 

fixed  z<0  the  decay  is,  however  exponential. 

(ii)  §«cond_order_solution__(  the_  f  ree_wayes_gotential| 

To  estimate  the  order  of  $2  we  again  retain,  like  in 

u 

Section  II2(ii),  only  w^  in  the  expression  of  w^  (2.23). 

For  f  given  by  (2.31)  we  immediately  obtain 


u 

W1 


r  0  _  0  a  1 

•r—  /  dx  /  dz  . . - . . 

**  -1  -h  ((x-x)2  +y2  +  (z+z)2] 


1  [  (x2+y2  ) 1//2  -  x)  {  [  (x+l)2  +  y2+h2 ) 1//2  -  x-1} 

"  57  <2-39) 

Based  on  the  results  of  the  preceding  section  we  can  retain  in 
(2.39)  only  the  term  originating  from  the  singularity  at  x»0, 
for  an  infinite  draft,  in  order  to  estimate  the  term  of  lowest 
order  in  F  ,  i.e . , 


'V'  yjj-  In  ( (x2+y2 ) 1,/2  -  x] 


(2.40) 


P2  and  4>2  have  the  expressions  (2.12)  and  (2.10),  respec¬ 


tively.  The  FT  of  w^(3w^/3x)  has  the  following  estimate 


(x  +yz ) 


/  dX  /  (In  r(l-cosX) JeirpcOB (X"6)/F*dr  (2, 


-Ti  o 
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where  the  last  expression  in  (2.41)  has  been  obtained  by  sub 
stituting  polar  coordinates  in  the  x,y  and  a,  6  planes.  The 
order  of  magnitude  of  (2.41)  may  be  estimated  from  the  In  r 
term 


(w 


uK 
i  vr 


) 


'V/ 


1 

57^ 


/  I 


-V  o 


(In  r  ) 


eir  p  cos (X-9)/F2  dr 


1  F2 

W  p 


(C  +  In  r 


+  In  |^r) 


(2.42) 


where  C  is  the  Euler  constant  (Gradshteyn  &  Rhyzhik,  1965). 
By  (2.42)  and  (2.10)  we  have 


♦2  % 


F6  (C  +  In  r  +  In  p-  2  In  F  ) 
pz (p-scc* S  -  iu  sec6) 


sec2 ft  e 


PZ/F2 


(2.4  3) 


Finally,  the  estimate  of  the  lowest  order  term  of  <t>~  becomes 
(2.5) 


$2  ^  F2ln  F  Im  /  d6 
1  -ir 


£-i(x  sec  0 +y  sij»  6  sec2 0+iz  sec28)/F2  ^  44) 


(lii)  DiBoueeion  o£  results 

Inspection  of  e4>^  (2.38)  and  e2^  (2.44)  shows  that  the 

amplitude  functions  are  0(eF2)  and  0(e2F2ln  F) ,  respectively. 
Hence,  the  expansion  of  the  free  waves  potential  is  not  uniform 
for  fixed  e  and  F-*-0,  although  the  singularity  is  much  milder 
than  that  corresponding  to  a  blunt  submerged  body.  To  ensure  the 
validity  of  the  thin  body  solution  the  condition  |e  In  F|  =  o(l) 
must  be  satisfied,  i.e.,  | B’ In [U ’/ (gL’ ) ]/L‘ |  <<  1,  where  L* 

is  the  forebody  length  (it  is  worthwhile  to  mention  that  this 
criterion  is  net  different  of  that  valid  for  two-dimensional  flow 
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past  a  body  of  similar  shape,  Table  1).  This  criterion  is  only 
marginally  satisfied  by  typical  slow  commercial  ships.  It  is 
emphasized,  however,  that  this  is  an  asymptotic  estimate  and  the 
actual  ratio  between  the  waves  amplitude  may  be  evaluated  only 
from  detailed  computations  of  and  •  Obviously,  the  above 

condition  applies  also  to  the  expansion  of  the  coefficient  of  wave 
resistance  (2.16). 

Again,  similarly  to  the  two-dimensional  flow,  blunter  bow 
shapes  will  impose  more  stringent  criteria  of  validity  of  the  thin 
body  solution  for  small  F  . 


4.  Derivation  of  Uniform  Small  Froude  Number  Solutions 


(i)  Velocity  straining 

The  argument  is  similar  to  that  given  in  Section  15:  we 
assume  that  the  presence  of  the  body  causes  a  free-surface  velocity 
straining  which  has  to  be  incorporated  in  the  first  order  solution 
since  it  causes  a  change  of  the  wave  number  and  not  only  of  the 
wave  amplitude. 


If  we  assume  that  all  the  terms  of  P2  (2.3),  (2.11)  contri¬ 
bute  to  the  F  lowest  order  terms  of  <J>2  ,  the  simple  straining 
of  the  horizontal  uniform  velocity  of  Section  15,  is  not  sufficient. 
The  generalized  straining  suggested  by  the  free  surface  conditions 
yields 


f2(1  +  3cu1>  f!c(2v1  !^y  +  “l  F!“l  IAt1 


(z  *  0,  y  0) 


(2.45) 
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(2.45)  is  supposed  to  replace  both  (2.2)  and  (2.3);  we  are  going 
to  prove  that  only  for  c  *=  o(l)  and  F  fix*d  it  does  indeed 
separate  into  these  two  equations. 

Keeping  the  body  boundary  condition  in  its  first  order  ver¬ 
sion#  we  write#  like  in  (2.7) 

j.  _  a (  Au  , 

<P  *  +  <P 

where  4>r  is  harmonic  and  regular  for  z<Q.  Substituting  (2.7) 
into  (2.45)  and  using  the  FT  relationships  (2.4)  and  (2.6)  we 
obtain  from  (2.45) 


$r(a,8;F) 


oo  ao 

TrfTr^-p-Iuar  /  dv  /  dT  K(a#8#v#x;F)  $r(v,T;F)  =• 

— oo  -OQ 


2F2w^(a,8) 

St  ■  .■ 

a^-p-iua 


(2.46) 


where 


K  =  (3(oc-v)a2  +  2(8-x)a8  -  2(a-v)za3  -  (ot-u)a‘‘]^(a-v#8-T)  .  ,  (2.47) 


(a# 8?F) | 2=0  satisfies#  therefore#  the  Fredholm  integral 
equation  (2.40).  This  equation  degenerates  precisely  into  the 
FT  of  the  thin  body  free  surface  conditions  (2.2)  and  (2.3)  if 
the  integral  equation  is  solved  by  successive  approximations. 

This  is  legitimate,  however,  only  for  sufficiently  small  values 
of  the  combination  between  e  and  F  which  multiplies  the 
second  term  of  (2.46);  this  combination  depends  on  the  order 
of  K  ,  which  in  turn  depends  on  the  order  of  .  A  sufficient 
condition  is,  however,  c  =  o(l)  F  =  0(1). 
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In  those  cases  in  which  F  is  so  small  that  the  usual  thin 
body  expansion  is  not  uniform,  $r  ,  solution  of  (2.46)  is  pre¬ 
sumably  a  valid  approximation  of  the  exact  solution  for  small 
and  F  . 

Solving  (2.46)  is  an  extremely  difficult  task.  We  consider, 
therefore,  the  method  of  coordinate  straining  as  an  alternative 
simplified  approach. 

( it )  Coordinate training 

Again,  the  argument  is  the  same  as  in  the  case  of  two- 
dimensional  flow:  it  is  assumed  that  the  nonlinear  free  surface 
effect  manifests  in  a  straining  of  the  coordinates  which  results 
in  a  virtual  displacement  of  the  system  of  singularities  represen¬ 
ting  the  image  of  the  body. 

We  begin  the  discussion  of  the  method  wit)  the  case  of  thin 
ships.  The  first  order  free  waves  potential  may  be  written  as 
follows  (2.22),  (2.35) 

*  *  tt/2 

*  Im  { jf  dxdz  /  exp{i( (x-x)sec  0  +  y  sin  0  eec20  - 

1  S  -tt/2 

-  i  (z+z)secJ0)/Fa  }sec20  d0)  .  (2.48) 


Like  in  (1.49)  we  may  assume  now  that  in  the  exponential 
function  of  (2.48)  x  and  z  are  replaced  by  x+  6x  and  z  +  6z, 
respectively,  where  6x,  <5z  are  straining  functions  of  order  e 
which  depend  on  x,  z  and  F,  Such  a  general  and  complicated 
straining  is,  nowever,  unnecessary;  fxora  the  discussion  of  Sec. 

113  it  is  seen  that  in  order  to  ensure  uniformity  for  F  -►  0  it 
is  enough  to  consider  the  displacement 
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of  the  moat  singular  points  of  3f(x,0)/3x  ,  x.v*.#  in  the  case  of 

finite  jumps  of  3f/3x  ,  the  points  of  abrupt  angle  changes  at 
the  water-line.  Obviously#  for  z«0#  6z  -  0  and  the  straining 
is  horizontal  solely.  Along  these  lines  let  us  consider#  for 
instance#  the  singularity  associated  with  the  bow  for  a  cylindri¬ 
cal  ship  of  finite  draft.  First#  before  straining,  <t>^  becomes 
by  integration  by  parts 


♦u>"  I 


JiF*  8f(V0)  ,*/2 


3x  -tt/2 


/  expfi  [  (x-x^sec  8 +y  sin  6  sec2 0- 


-iz  sec^l/F2  }cos  0  du 


(2.49) 


where  x^  is  the  abscissa  of  the  bow  for  z*0;  for  the  sake  of 
simplicity  we  may  take  (like  in  Sec.  113)  3f(xt#0)/3x  »  -1  which 

is  tantamount  to  defining  e  as  the  tangent  of  the  entrance  angle. 
Assuming  now  that  the  free  surface  nonlinear  effect  leads  to  a 
straining  6*b  '  ♦lb  becomes 


w  21F2  ^ 

^lb  “  ~Irn  “ —  /  exp(i  [  (x-x^-ix^)  sec  0  +y  sin  0  sec20  - 
-n/2 

-  iz  sec28]/F2}  oos  0 d0 


(2.50) 


Like  in  Sec.  15  ix^  =  0(e)  is  sought  by  requiring  that  for 
fixed  F  and  e  =  o(l),  the  thin  body  expansion  «  e<>^H  +  c 2  <J>^ 

of  the  free  waves  potential  remains  uniform  as  F-*-0,  i.e.#  the 

second  order  approximation  is  not  more  singular  them  the  first. 


The  expansion  of  (2.50)  for  (x-x^)2  +  z2  ^  0  yields 


.w  _.w  _  24xb  rn 


^  /  exp(i[  (x-x^sec  0  +  y  sin  0  sec20  -  iz  sec^J/F^dS 

“n/2  (2.51) 


HYDRONAUTICS,  Incorporated 

43. 

6xb  has  now  to  be  determined  such  that  the  .last  term  of 
(2.51)  should  cancel  the  lowest  F  order  term  of  4)^  .  It 
turns  out  that  this  is  possible  for  the  estimate  of  4^  given 
in  (2.44).  As  result  Sx^  =  0(eFJln  F)  which  is  the  same  as 

in  the  two-dimensional  case  (Table  1).  Sx^  can  be  substituted 

in  which  becomes  a  valid  approximation  even  if 

cln  F  =  0(1).  Obviously,  the  nonuni formity  of  the  thin  body  usual 
approximation  in  the  latter  case  is  a  result  of  the  expansion  of 
the  wave  number  function  in  (2.50)  in  an  c  power  series. 

The  expression  of.  the  tree  waves  (2,50)  differs  from 

4>lb  (2.49)  only  by  a  change  cf  phase  resulting  from  an  apparent 
displacement  of  the  bow,  which  i.s  velocity  dependent.  The  same 
is  true  for  the  curve  representing  the  coefficient  of  wave  resis¬ 
tance  as  function  of  F  .  This  result  is,  at  least  in  principle, 
in  agreement  with  experimental  finding- . 

If  there  are  other  points  of  slops  discontinuity,  additional 
^training  factors  have  to  be  incorporated  in  ^  »  If  the  pro¬ 
cedure  is  applied  to  a  submerged  body  there  is  ar.  additional  free¬ 
dom  in  selecting  a  vortical  straining;  the  choice  between  horizon¬ 
tal  Kc\d  vertical  strainings  depends  on  whether  we  have  to  cancel  P 
or  Q  type  functions  in  the  expression  of  R  (2.14).  In  the 
cane  of  a  source-like  blunt  nose  (Sec.  112)  the  straining  is  of 
order  c/F  as  compared  to  the  order  e  in  two-dimensional  flews 
(Table  1) .  | 

Since  our  estimate  of  the  lowest  order  term  of  4^  (2.44)  } 

was  based  only  on  part  of  the  expression  of  (x,y)  (2.3),  it  j 

is  not  sure  that  the  complete  wave  spectrum  functions  P^(6)  l 

w  *  I 

and  of  4>2  are  the  same  as  in  (2,44).  If  they  are  different  | 

the  coordinate  straining  m&y  be  successful  only  if  we  assume  that 
6x^  is  a  function  of  b  ,  which  makes  the  straining  less  meaningful  S 

then  in  the  case  of  two-dim^ nsional  flow.  This  question  is,  how¬ 
ever,  iefv  open  at  present. 
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5.  Conclusions 

The  results  obtained  in  the  study  of  the  small  r'roude  numbo r 
limit  of  the  thin  body  expansion  in  the  case  of:  three-dimensional 
flows  do  not  differ  in  principle  from  those  pertaining  to  two- 
dimensional  flews;  the  computations ,  however,  ere  much  more  tedrou 
Again,  the  non uniform  behavior  of  the  expanse  on  of  the  potenti  al 
of  the  free  waves  is  apparently  related  to  the  crc-seno-.?  of  the 
small  parameter  e  in  the  wave  number  of  the  spectrum  function, 
and  not  only  in  the  amplitude.  it  is  worthwhile  to  mention  th.vu 
the  nonlinear  effects  considered  here  are  associated  with  tentc  of 
the  second  ord»r  pressure  which  result,  from  tie  local  disturbance 
of  the  free  -surface.  The  interaction,  between  the  first  and  second 
order  free  waves  yielda  terms  which  are  or:  higher  order  as  F-0. 

A  relatively  simple  method  oi  rendering  uniform  the  expansion 
of  the  free  waver  potential  suggested  here  is  the  coordinate 
straining.  It  result,’;,  ir.  an  apparent  horizontal  displacement  of 
the  singularities  of  the  water-line  contour  in  the  case  of  thin 
ships.  This  displacement  is  of  the  order  of  ti»e  beam  and  its 
dependance  on  the  Freuds  number  is  related  to  th  -  nature  of  the 
singularity.  The  actual  value  of  the  straining  factor?  has?  to 
be  computed  numerically  for  each  particular  case.  The  influence 
of  the  straining  becomes  appreciable  tor  blunt  bodice  moving  at 
low  speeds. 
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The  expression  of  f2  (1.21)  beoomes,  by  substitution  of 
ux  (1.22)  and  vx  (1.23) 

tt) 

f2  i  /  do  w(^r){-y  IwJ(o)  -wj(c)]2  -  [wj(o)  -w1(o)]  * 

x  [  /  dx  wj'(x)  w(2-y-))  -  /  dx^  /  dx2  fw^(x^)w^(T2)u(-^jri)u(-p3r2-)  + 

— oo  ..co  »ao 

9  „  -  °‘T'  °~Ty  f  .  O-X  _  o-x 

+  2w,  (x.,  )w,  (x,,)  w  -  Hpt- )  -4  %>,  (x,)w,  (t5)  m (i^r  )  u («y«*i  )  }  (A.l) 

x  JL  *  a  a  x  1  1  ^  r  r 


By  integration  by  parts  it  can  be  shown  that 


00  cr-x  o-x-  »  .  a-x.  a-x 

/  a) (-gr^) « )  (-|rr)  do  “  IF2  /  — a) (■  p— ) )  <D(-£y)do  (A.  2) 

BO  "OO  1  2 

00  G~T  CT“T0  00  .  Q-T*  .  0—*[\ 

/  ut-pi)  (-p?^  <^r)  do  -  -iF2  /  3^-«(-T7i)+3iru(-1pi)]  u(^)do  (A.  3) 

60  — OO  1  2 


By  residues  we  also  have 
00  wY1  ( x0 ) 

/  -i — 2-  dx0  «  2ni  w?(o) 
'  U  -  X-  2  1 

-®  i 


r  wi(ti) 

/  -± - =—  dx 

L  o-x. 


»  -2Tri  w^o)  (A. 4) 


Substituting  (A, 2),  (A. 3)  and  (A. 4)  in  (A.l)  results  in  the  final 
expression  of  fj(z)  (1.25). 
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The  complete  expression  of  f 2  is  given  in  (1.29).  By 
residues  we  obtain  for  the  last  integral  of  (1.28) 


/  (— - - —“>  w(?^£>)  dp  =  2TTi  [u{LlL)  -uflli)] 

-oo  p“Z„  P"Z.  F2  F2  F2 

JC  t 


c-z. 


(A.  5) 


We  represent  now  u>  like  in  (1.18) 


°~zt  *iX'/F 

w( — -i)  «  /  - - 

F2  o  X  +  o  -  z, 


dA 


(A. 6) 


and  similarly  u(a-zt/F2).  To  derive  the  expression  of^the  far 
free  waves  we  replace  w(z-a/F2)  by  2ni  e-i2//F  ei<7//F  .  With 

these  transformations  the  lest  term  of  f2  becomes 


i  iA/F2  t  ,  1 


/ 


TTlP2  O 


/  (• 


i— )  eio/ri  do  (A.  7) 


-0°  °“z£  °~zfc  0"zt  A+a-z^  A+a-zt 


Again,  by  using  the  residue  theorem  in  the  last  integral,  (A. 7) 
becomes 


nr 


2  L^l^2  l 

V<o  1  /  n-T= 


1  1  .  1  ,  ,  e 


iA/F2 


+  — 3dx~ 


-  e 


-y*-2  »eiX/F 


/ 


7=T 


dA) 


(A. 8) 


^  i 


where  h  -  im  it  -  fc,  ^  and  n*  I(  ,  1(  R,  ^  ^ 

Integrating  in  (a. 8)  and  adding  the  reeiduea  of  the  first 
tent  of  f2U.28)  lead^  to  the  final  expression  of  f“  a. 32) 
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FIGURE  1  -  TWO-DIMENSIONAL  FLOW  PAST  A  SUBMERGED  BODY 


